Generalised constraints on the curvature perturbation from primordial black holes 
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Primordial black holes (PBHs) can form in the early Universe via the collapse of large density 
perturbations. There are tight constraints on the abundance of PBHs formed due to their gravita- 
tional effects and the consequences of their evaporation. These abundance constraints can be used 
to constrain the primordial power spectrum, and hence models of inflation, on scales far smaller than 
those probed by cosmological observations. We compile, and where relevant update, the constraints 
on the abundance of PBHs before calculating the constraints on the curvature perturbation, taking 
into account the growth of density perturbations prior to horizon entry. We consider two simple 
parameterizations of the curvature perturbation spectrum on the scale of interest: constant and 
power-law. The constraints from PBHs on the amplitude of the power spectrum are typically in the 
ranee 10~ 2 - 10" 1 with some scale dependence. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 



Primordial black holes (PBHs) can form in the early 
Universe via the collapse of large density perturba- 
tions [HQ- If the density perturbation at horizon en- 
try in a given region exceeds a threshold value, of or- 
der unity, then gravity overcomes pressure forces and 
the region collapses to form a PBH with mass of order 
the horizon mass. There are a number of limits, span- 
ning a wide range of masses, on the PBH abundance. 
PBHs with mass Mpbh ^$5 x 10 14 g will have evapo- 
rated by the present day H 0| and their abundance is 
constrained by the consequences of the Hawking radia- 
tion emitted. More massive PBHs are constrained by 
their present day gravitational effects. The resulting lim- 
its on the initial mass fraction of PBHs are very tight, 
(3 = ppBu/ptot < 0(1O~ 20 ), and can be used to con- 
strain the power spectrum of the primordial density, or 
curvature, perturbations (see e.g. Ref. @). 

The power spectrum of the primordial curvature per- 
turbation, Vn(k), on cosmological scales is now accu- 
rately measured by observations of the cosmic microwave 
background (CMB) Q and large scale structure @, [I]- 
These measurements can be used to constrain, and in 
some cases exclude, inflation models (c.f. Ref. [9]). Cos- 
mological observations span a relatively small range of 
scales (comoving wavenumbers between k ~ 1 Mpc -1 and 
k ~ 10 -3 Mpc -1 ), and hence probe a limited region of 
the inflaton potential. The PBH constraints on the cur- 
vature power spectrum are fairly weak; the upper limit 
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is many orders of magnitude larger than the measure- 
ments on cosmological scales. They do, however, apply 
over a very wide range of scales (from k ~ 10 _2 Mpc _1 
to k ~ 10 23 Mpc _1 ) and therefore provide a useful con- 
straint on models of inflation [l(| ■ 

The simplest assumption for the power spectrum is a 
scale-free power law with constant spectral index, n: 



Vn{k) 



2^ 



(\K k 



Vn(k ) 



(1) 



where k is a suitably chosen normalisation scale. In 
this case the PBH abundance constraints require n < 
1.25 - 1.30 [H El El 03 ■ The spectral index on cos- 
mological scales is, however, now accurately measured: 
n = 0.9631q pig 6]. In other words, if the power spec- 
trum is a pure power law then the number of PBHs 
formed will be completely negligible. However, as we will 
now outline, if the primordial perturbations are produced 
by inflation then the power spectrum is not expected to 
be an exact power law over all scales. 

The power spectrum produced by slow-roll inflation 
can be written as an expansion about a wave number kg 
(e.g. Ref. [H) 



lnV n (k) ra lnP K (feo) + [n(fc )-l]ln 



k 



2 \k 



(2) 



where the spectral index and its running, cv(fc) = 
dlnn/dlnfc, are evaluated at fco and can be expressed 
in terms of the slow-roll parameters. This expansion is 
valid provided In (fc/fcn) is small for the relevant k values. 
This is the case for cosmological observations, but not for 
the wide range of scales probed by PBH constraints. 
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In fact only very specific, and contrived, inflaton poten- 
tials produce a constant spectral index [16[ . It is possible 
for the power spectrum to vary sufficiently with scale so 
that PBHs can potentially be over-produced. For in- 
stance in the running-mass inflation model [l7l fl8| , the 
power spectrum is strongly scale-dependent and PBH 
constraints exclude otherwise viable regions of parame- 
ter space pjj ■ More generally, Peiris and Easther [2l[ 
recently found, using slow roll reconstruction, that infla- 
tion models which are consistent with cosmological data 
can over-produce PBHs. 

Motivated by this, we compile and update the con- 
straints on the abundance of PBHs (Sec. [TTJ) . We then 
translate the abundance constraints into detailed gener- 
alised constraints on the power spectrum of the curvature 
perturbations (Sec. IHIj) . taking into account the evolu- 
tion of the density perturbations prior to horizon entry 



II. PBH ABUNDANCE CONSTRAINTS 

The PBH constraints can, broadly, be split into two 
classes: those that arise from their present day gravita- 
tional consequences and those that arise from the prod- 
ucts of their evaporation. In both cases, in order to con- 
strain the primordial density or curvature perturbation, 
we need to translate the constraints into limits on the 
initial PBH mass fraction. 

Throughout we will assume that the PBHs form at a 
single epoch and their mass is a fixed fraction, /m , of the 
horizon mass M PB h = ImM h , where f M ~ (1/3) 3 / 2 [HI]. 
A scale invariant power spectrum produces an extended 



PBH mass function dnpsH/d-MpBH oc M, 



-5/2 
PBH 



@, m, 



however (as discussed in the introduction above) in this 
case the number density of PBHs would be completely 
negligible [l(| For scale-dependent power spectra 

which produce an interesting PBH abundance it can 
be assumed that all PBHs form at a single epoch [25| . 
As a consequence of near critical phenomena in gravi- 
tational collapse [26|, [27], Hl[ the PBH mass may, how- 
ever, depend on the size of the fluctuation from which it 
forms l3lj in which case the mass function has 

finite width. Most of the constraints that we discuss 
below effectively apply to the mass function integrated 
over a range of masses. The range of applicability is usu- 
ally significantly larger than the width of the mass func- 
tion produced by critical collapse, so in the absence of 
a concrete prediction or model for the primordial power 
spectrum in most cases it is reasonable to approximate 
the mass function as a delta-function. The constraints 
from cosmic-rays and gamma-rays produced by recently 
evaporating PBHs are an exception to this. These con- 
straints depend significantly on the PBH mass function 
and therefore need to be calculated on a case by case ba- 
sis [H], [33|, H3, [H, [36[ . We therefore do not include these 
constraints in our calculation of generalised constraints 
on the curvature perturbation power spectrum. 

Taking into account the cosmological expansion, the 



initial PBH mass fraction, /3(Mpbh), is related to the 
present day PBH density, ^ PBH , by 
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(3) 

where a is the scale factor, 'eq' refers to matter-radiation 
equality and p cr it is the critical energy density. Using 
the constancy of the entropy, (s = g* s a 3 T 3 ), where g± s 
refers to the number of entropy degrees of freedom, we 
relate the scale factor to the temperature of the Universe 
and using the radiation density, p — ^g+T , and horizon 
mass, Mh = ^-pH~ 3 , we obtain 
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(4) 



where g+ is the total number of effectively massless de- 
grees of freedom and we have taken g+ s ss <?*. The hori- 
zon mass at matter-radiation equality is given by (c.f. 
Ref. [13) 



M oq -—p H- 3 - 



P°ad 
3 G-eq &cq 



(5) 



Inserting numerical values, g\ ~ 100, gl q ~ 3, Q® ad h 2 
4.17 x 10~ 5 , 
0.07 n° m h 2 Mpc" 1 , 



0.1326 ± 0.0063 
so that 



1.88 x lO-^gcm-^, fc oq = 
= 24000 fl° m h 2 and fl^h 2 = 
gives = 1.3 x 10 49 (fl m h 2 )- 2 g 



/3(M PBH ) = 6.4 x 10~ 19 n' 



PBH 
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As first shown by Hawking [38] a black hole with mass 
Mpbh emits thermal radiation with temperature 



1 



/10 13 g\ 



w 1.06 GeV 

SttGMpbh \ Mpbh J 



(7) 



The current understanding of PBH evaporation [39j is 
that PBHs directly emit all particles which appear ele- 
mentary at the energy scale of the PBH and have rest 
mass less than the black hole temperature. Thus if the 
black hole temperature exceeds the QCD confinement 
scale, quark and gluon jets are emitted directly. The 
quark and gluon jets then fragment and decay produc- 
ing astrophysically stable particles: photons, neutrinos, 
electrons, protons and their anti-particles. Taking into 
account the number of emitted species the mass loss rate 
can be written as 14011 



dM 



PBH 



dt 



-5.34 x 1O 25 0(M] 



PBH 



i^pbhS 8 1 



(8) 



where 0(Mpbh) takes into account the number of di- 
rectly emitted species (4>(Mpbh) — 0.267go + 0.147gi/ 2 + 
0.06gi + 0.02^3/2 + O.OO752 where g s is the number of 
degrees of freedom with spin s) and is normalized to 
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one for PBHs with mass Mpbh ^ 10 g which can 
only emit photons and neutrinos. For lighter PBHs 
0(5 x 10 14 g < Mpbh < 10 17 g) = 1.569. The PBH life- 
time is then given by [40| 

r » 6.24 x 1O- 27 M| bh 0(A/pbh) _1 s . (9) 

From the WMAP 5 year data the present age of 
the Universe is to = 13.69 ± 0.13 Gyr. The initial 
mass of a PBHs which is evaporating today is therefore 
Mpbh ~ 5 x 10 14 g jij, while less massive PBHs will have 
evaporated by the present day. 

We will now compile, and where relevant update, the 
PBH abundance constraints. We divide the constraints 
into two classes: those, for PBHs with Mpbh > 5 x 10 14 g, 
arising from their gravitational consequences (Sec. Ill A|) 
and those for Mpbh < 5 x 10 14 g arising from their evap- 
oration (Sec. lHB) . 

A. Gravitational constraints 



tern [4^|. A null search using BATSE data leads to a 
constraint [44| 

Q c < 0.2 for 1(T 16 M < Mpbh < 1O~ 13 M , (12) 

where fi c is the density of compact objects, assuming a 
mean GRB red-shift of one. 



Quasars 

Compact objects with mass 10~ 3 M Q < Mpbh < 
300Af Q can microlens quasars, amplifying the continuum 
emission without significantly changing the line emis- 
sion . Limits on an increase in the number of small 
equivalent width quasars with red-shift lead to the con- 
straint 

fl c < 0.2 for O.OOIA/q < M PBH < 60M Q , (13) 
assuming il tot = fl c - 



1. Present day density 

The present day density of PBHs with Mpbh > 5 x 
10 14 g which haven't evaporated by today must be less 
than the upper limit on the present day cold dark mat- 
ter (CDM) density. Using the 5 year WMAP measure- 
ments d, ft^rw/i 2 = 0.1099 ± 0.0062, h 



71O+0.026 
u - ' lt, -0.027' 



gives 



which, using Eq.©, leads 



ftp BH < 0.25 . 



(10) 



to 



f3(M } 



W<i..xio- 

for Mpbh > 5 x 10 14 g . 



2. Lensing of cosmological sources 



1/2 



(11) 



If there is a cosmologically significant density of com- 
pact objects then the probability that a distant point 
source is lensed is high ,41!]. The limits as given below 
have been calculated assuming an Einstein de Sitter Uni- 
verse, fi m = 1, and a uniform density of compact objects. 
The recalculation of the constraints for a A dominated 
Universe would be non-trivial. The constraints would, 
however, be tighter (due to the increased path length 
and the larger optical depth to a given red-shift) [42l |. 
and the constraints given below are therefore conserva- 
tive and valid to within a factor of order unity. 



Radio sources 

Massive compact objects, 10 6 M Q < Af PB H < 10 8 M Q , 
can millilens radios sources producing multiple sources 
with milliarcsec separation [46j . A null search using 300 
compact radio sources places a constraint |47| 



fl c < 0.013 for 1O 6 M < Mpbh < 10 8 M G 



3. Halo fraction constraints 



(14) 



There are also constraints from the gravitational con- 
sequences of PBHs within the Milky Way halo. They 
are typically expressed in terms of the fraction of the 
mass of the Milky Way halo in compact objects, /h = 
Mp^/M^w. They require some modeling of the Milky 
Way halo (typically the density and/or velocity distribu- 
tion of the halo objects). Consequently there is a factor of 
a few uncertainty in the precise values of the constraints. 

Assuming that PBHs make up the same fraction of 
the halo dark matter as they do of the cosmological cold 
dark matter, and ignoring the uncertainties in the CDM 
density (since this is negligible compared with the uncer- 
tainties in halo fraction limit calculations), we can relate 
the halo fraction to the PBH cosmological density: 



fh 



M 



MW 
PBH 



M MW 
M CDM 



Ppbh 

PCDM 



"PBH ft 



5Ql 



PBH 



(15) 



Microlensinq 

Gamma-ray bursts 

Solar and planetary mass compact objects in the Milky 
Light compact objects can femtolens gamma-ray bursts Way halo can microlens stars in the Magellanic Clouds, 
(GRBs), producing a characteristic interference pat- causing temporary one-off brightening of the microlensed 



4 



star [48|. The relationship between the observed optical 
depth, t, (the probability that a given star is amplified 
by more than a factor of 1.34) and the fraction of the halo 
in MACHOs depends on the distribution of MACHOS in 
the MW halo. For the 'standard' halo model used by 
the microlensing community (a spherical cored isother- 
mal sphere) r«5x 10~ 7 /h [4?J with the derived 
value of limits on fh varying by factors of order unity for 
other halo models. 

The EROS collaboration find a 95% upper confidence 
limit t < 0.36 x 10~ 7 which they translate into limits on 
the halo fraction [5l| : 

f h < 0.04 for 10- 3 M Q < Mpbh < 1O _X M , (16) 



disk stars [58|. This leads to a limit, from the observed 
stellar velocity dispersions, on the halo fraction in mas- 
sive objects [59j 



fh < 



Mm 



disk.lim 



A/disk, lim 
ATpBH 

3 X 10 6 

Cob 



Ph 



0.01M Q pc- 3 

2 



60kms 



1Q1° yr 



M , (25) 



where ph is the local halo density and er b s and t s are the 
velocity dispersion and age of the halo stars, respectively. 



f h < 0.1 for 10- 6 M Q < A/pbh < M { 



- 



(17) 



Combined EROS and MACHO collaboration limits on 
short duration events constrain the abundance of light 
MACHOs (H 



f h < 0.25 for IQ- 7 M Q < M PBH < 10~ 3 M ( 



(18) 



while a dedicated search by the MACHO collaboration 
for long (> 150 days) duration events leads to limits on 
more massive MACHOs [H: 



fh < 1.0 for O.3M < A/pbh < 3OM , (19) 

f h < 0.4 for a/pbh < 10M Q . (20) 
Combined, these limits give 



or 



f h < 0.25 for 10" 7 A/ Q < M PBH < 10~ 6 Af Q , (21) 
f h < 0.1 for 10- 6 Af o < A/pbh < M Q , (22) 
f h < 0.4 for Af < A/pbh < 1OM . (23) 

Wide binary disruption 

More massive compact objects would affect the orbital 
parameters of wide binaries [54], [55| . Comparison of the 
separations of observed halo binaries with simulations of 
encounters between compact objects and wide binaries 
lead to a constraint [56f 



f h < 0.2 for 10 3 Af o < A/pbh < 10 s M, 







(24) 



See however Ref. [57J for a recent re-examination of this 
constraint which leads to a somewhat weaker limit. 



Disk heating 

Massive halo objects traversing the Galactic disk will 
heat the disk, increasing the velocity dispersion of the 



B. Evaporation constraints 

1. Diffuse gamma-ray background 

PBHs with masses in the range 2 x 10 13 < M P bh < 5 x 
10 14 g evaporate between z w 700 and the present day and 
can contribute to the diffuse gamma-ray background [H, 
011,011 \&\ . As discussed above, these constraints 
depend significantly on the PBH mass function and hence 
we will not consider them further. 

2. Cosmic-rays 

The abundance of PBHs evaporating around the 
present day can also be constrained by limits on the 
abundance of cosmic-rays (in particular positrons and 
antiprotons) (23l . [63| . The constraints from anti-protons 
have been calculated for several mass functions and are 
essentially equivalent to those from the diffuse gamma- 
ray background [H, [64| . 

3. Neutrinos 

Neutrinos produced by PBH evaporation contribute to 
the diffuse neutrino background. The neutrino spectrum, 
and hence the resulting PBH abundance constraints, de- 
pend strongly on the PBH mass function, but the con- 
straints are typically weaker than those from the diffuse 
gamma-ray background [33j, l34j . 

4- Hadron injection 

PBHs with mass Mpbh < 10 10 g have a lifetime 
t < 10 3 s and evaporate before the end of nucleosyn- 
thesis, and can therefore affect the light element abun- 
dances [H, [gl, [6?], The constraints from hadron 
injection have been re-evaluated (see Ref. [6!|), taking 
into account the emission of fundamental particles and 
using more up to date measurements of the Deuterium 
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and 4 He abundances (D/H < 4.0 x 10~ 5 ,Y P < 0.252 
respectively) : 

P(Mp BH ) < 10- 20 for 10 8 g < Mpbh < 10 10 g,(26) 
/3(M P bh) < HT 22 

for 10 10 g < A/pbh < 3 x 10 10 g. (27) 

5. Photodissociation of deuterium 

The photons produced by PBHs which evaporate be- 
tween the end of nucleosynthesis and recombination can 
photodissociate deuterium [7C| . The resulting constraints 
on the PBH abundance have been u pda ted, in the con- 
text of braneworld cosmology in Ref . [7l| ■ Adapting that 
calculation to the standard cosmology we find: 

Ipbh < 10 13 g. (28) 



for 10 10 g < M P 



6. CMB distortion 

Photons emitted by PBHs which evaporate between 
z ~ 10 6 and recombination at z ~ 10 3 can produce 
distortions in the cosmic microwave background radia- 
72|. Using the COBE/FIRAS limits on spectral 



tion 

distortions of the CMB from a black body spectrum [73 
Ref. [zi finds 

/3(M PBH ) < KT 31 for 10 11 g < M PBH < 10 13 g . (29) 

7. ( Quasi- )stable massive particles 

In extensions of the standard model there are gener- 
ically stable or long lived massive (0(100 GeV)) parti- 
cles. Light PBHs with mass Mpbh ^ 10 11 g can emit 
these particles and their abundance is hence limited by 
the present day abundance of stable massive particles [75j ] 
and the decay of long-lived particles [7rl l77l] 1 . 

Gravitinos in supergravity theories and moduli in 
string theories are generically quasi-stable and decay af- 
ter big bang nucleosynthesis, potentially altering the light 
element abundances. The effect of their decay on the 
products of big bang nucleosynthesis leads to a constraint 
on the initial PBH fraction [76j |: 



/3(M PBH ) < 5 x 10- 19 (A 



1/4 



(?) 



v 6 x 10" 3 / 



/mM] 



PBH 



-1/2 



10 9 g J \1Q- 14 
for Mpbh < 10 9 g, 



(30) 



where is the fraction of the luminosity going into 
quasi-stable massive particles, g\ is the initial number of 
degrees of freedom (taking into account supersymmetric 
particles), a is the mean energy of the particles emit- 
ted in units of the PBH temperature and is the limit 
on the quasi-stable massive particle number density to 
entropy density ratio. 

In supersymmetric models, in order to avoid the decay 
of the proton, there is often a conserved quantum num- 
ber R-parity, which renders the Lightest Supersymmetric 
Particle (LSP) stable and the present day density of such 
stable particles produced via PBH evaporation must not 
exceed the upper limit on the present day CDM den- 
sity [75|]. This leads to a constraint on the initial PBH 
fraction (c.f. Ref. |z|): 

«M PB „)< 6 xio-^(A)" 4 (|)(: 



PBH 



10 



11 , 



200 

1/2 



^LSP 

0.6 



/ wlsp 
V 100 GeV 



for Mpbh < 10 11 g 



100 GeV\ 



(31) 



where ttilsp is the mass of the LSP and slsp is the frac- 
tion of the luminosity carried away by the LSP. 

These constraints depend on the (uncertain) details 
of physics beyond the standard model, and we therefore 
summarise them conservatively as 

/mMpbh \ 1//2 
10 n g ) 
for Mpbh < 10 11 g . (32) 



/3(M PBH ) < 10- 



8. Present day relic density 

It has been argued that black hole evaporation could 
leave a stable Planck mass relic @, [7^, H3|, in which 
case the present day density of relics must not exceed 
the upper limit on the CDM density 

n° cl < 0.25 . (33) 

Writing the relic mass as M re i = / ro iMpi this gives a 
tentative constraint 

1 / A/f — ^ 3 ^ 2 



PBH 



f 1 /^ f , \ 5 x 10 14 1 
for Mpbh < 5 x 10 14 g. (34) 



1 More massive PBHs can also emit these particles in the late 
stages of their evaporation, when their mass drops below ~ 10 9 g. 
However the resulting constraints are substantially weaker than 
those from hadron injection during nucleosynthesis. 



The constraints are summarised in table U and are dis- 
played in Fig. 1. As can be seen from Fig. 1, the con- 
straints probe a very large range of scales and in some 
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cases several constraints overlap across particular mass 
ranges. The solid line indicates the strongest constraints 
for each mass scale and we consider only these when con- 
straining the primordial power spectrum in Sec. IIIII 



III. CONSTRAINTS ON THE CURVATURE 
PERTURBATION POWER SPECTRUM 

We focus in the following on the standard case of PBH 
formation, which applies to scales which have left the 
horizon at the end of inflation. It has recently been 
shown [8l], H3] that PBHs can also form on scales which 
never leave the horizon during inflation, and therefore 
never become classical. We also only consider gaussian 
perturbations and a trivial initial radial density profile, 
and refer to Ref. [83| for the effects of non-gaussian per- 
turbations and to Refs. [UHll for estimates on the effect 
of deviations from a trivial initial density profile. 

A region will collapse to form a PBH if the smoothed 
density contrast, in the comoving gauge, at horizon 
crossing (R = (aH)~ 1 ), 5h m (R), satisfies the condition 
<5c < ^hor(-R) < 1 [24I) where S c ~ 1/3. The mass of the 
PBH formed is approximately equal to the horizon mass 
at horizon entry, Mpbh = JmMh, which is related to the 
smoothing scale, R, by [37j 



M H = M^(k cq R) 



9*, eg 

9* 



1/3 



(35) 



where M^ q = 1.3 x 10 49 (£l m /i 2 )~ 2 g is the horizon mass 
at matter-radiation equality. 

Taking the initial perturbations to be Gaussian, the 
probability distribution of the smoothed density contrast, 
P{6 hoT {R)), is given by (e.g. Ref. ||) 



exp 



^a hor (R) 
where cr(R) is the mass variance 



a 2 (R)= / W 2 (kR)V s (k,t) — , 



(36) 



(37) 



and W(kR) is the Fourier transform of the window 
function used to smooth the density contrast. We as- 
sume a Gaussian window function for which W(kR) = 
exp(-fc 2 i? 2 /2). 

The fraction of the energy density of the Universe con- 
tained in regions dense enough to form PBHs is then 
given, as in Press-Schechter theory [§7j by, 



/3(M PBH ) 



MpBH 

' M H 



Sc 



P(S hor (R)) d5 hm (R) . (38) 



This leads to a relationship between the PBH initial mass 
fraction and the mass variance, 



/3(Mpbh) = 



2/: 



M 




20 30 
logio M 

FIG. 1: The limits on the initial mass fraction of PBHs as a 
function of PBH mass (in grams). The solid lines represent 
the tightest limits for each mass range and the dotted lines 
are the weaker limits where there is an overlap between con- 
straints. As discussed in Sec. [TT] we have not considered the 
diffuse gamma-ray background constraint which applies for 
2 x 10 13 g < Mpbh < 5 x 10 14 g as it depends significantly on 
the PBH mass function. 



x / exp 

'<5c 



/ M erfc 



SI 



(R) 



Hor(R) 

Sc. 



ddhov(R) 



V2a hol (R) 



(39) 



The constraints on the PBH initial mass fraction can 
therefore be translated into constraints on the mass vari- 
ance by simply inverting this expression. 

In order to calculate the mass variance we need the 
density contrast in the comoving, or total matter, gauge 
as a function of time and scale (c.f. Ref. [3, [36[). To 
calculate this we take the expressions for the evolution of 
perturbations in the conformal Newtonian gauge, valid 
on both sub- and super-horizon scales, and carry out a 
gauge transformation to the total matter gauge (for fur- 
ther details see Appendix [A")) . We find 



ji(k/V3aH)Tl, (40) 



27rcr hor (i?) 



where j\ is a spherical Bessel function and 1Z is the pri- 
mordial curvature perturbation. Hence the power spec- 
trum of the density contrast is given by: 

n(M) = y(^) jf(k/V3aH)V K (k). (41) 
Substituting this into Eq. ([3T)) . and setting R = (aff) -1 , 
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TABLE I: Summary of constraints on the initial PBH abundance, /3(Mpbh)- 



description 


mass range 


constraint on /3(Mpbh) 


Gravitational constraints 


present day PBH density 


Mpbh > 5 x 10 14 g 


<2xlO- 19 ( 7I J^k i ) 


GRB femtolensing 


1O- 16 M < Mpbh < 1(T 13 M 


/ \ 1/2 
s i in— 19 1 Mpbh 1 

<1XlU {fMSxW^z) 


Quasar microlensing 


O.OO1M < AfpBH < 60M Q 


<lxl0- 19 ( 7 ^k i ) 


Radio source microlensing 


10 6 M Q < Mpbh < 10 s M Q 


<6xl0- 20 ( ?; ^k i ) 1/2 


Halo density a 


LMC Microlensing 


KT 7 M < Mpbh < 1O" 6 M 


<3xl0- 20 ( 7 ^&i) 

/ \ 1 / 2 

< 1 X 10 _2 ° ( , f F f,U4 

< 5 x 1(T 20 (-, \ 




KT 6 M < Mpbh < M 




M < Mpbh < 1OM 


Wide binary disruption 


1O 3 M < Mpbh < 10 s M 


f \ 1 l 2 
< 3xl °- 20 (77&i) 


Disk heating 


Mpbh > 3 x 1O 6 M 


<2xl0 6 -^(^f i )" /J 


Evaporation 


diffuse gamma-ray background 


2 x 10 13 g < Mpbh < 5 x 10 14 g 


depends on PBH mass function 


cosmic-rays 


similar to DGRB 


depends on PBH mass function 


neutrinos 


similar to DGRB 


depends on PBH mass function 


hadron injection 


10 s g < Mpbh < 10 10 g 
10 10 g < Mpbh < 3 x 10 10 g 


< 1(T 20 

< 1(T 22 


photodissociation of deuterium 


10 10 g < Mpbh < 10 13 g 


< 3 >< 10 ~ 22 


CMB distortion 


10 11 g < Mpbh < 10 13 g 


< 10~ 21 


(Quasi-)stable massive particles b 


Mpbh < 10 11 g 


/ \ -I/ 2 

^. . in-18 / f hi ^pbh 1 
<~ 10 1 lO", ) 


present day relic density c 


Mpbh < 5 x 10 14 g 


/ \ 3/2 
^ A 1 ( »fpBH 1 



"These constraints depend on the PBH distribution within the 
Milky Way halo and hence have a factor of order a few uncertainty. 

b Conservative summary, depends on physics beyond the standard 
model of particle physics. 

c Only applies if evaporation leaves stable relic. 



gives 

<r(R) = y J o (kRfjUkR/Vs) 

xexp(-fc 2 i? 2 )^(fc)^. (42) 

Since the integral is dominated by scales k ~ l/R we 
assume that, over the scales probed by a specific PBH 
abundance constraint, the curvature power spectrum can 
be written as a power law 

/ u \ n(*o)-l 

Vn(k)=Vn(k ){^\ . (43) 

This assumption is valid for general slow-roll inflation 
models such as those considered in Refs. pjl, [2(J l2lj . 



Using Eqs. ((35J) and (g2J| we can translate the PBH 
abundance constraints in Sec. [IT] into constraints on the 
amplitude of the curvature perturbation spectrum. For 
each constraint we take the pivot point, fcrjj to corre- 
spond to the scale of interest, k = l/R, and consider 
a range of values for n(fco) consistent with slow roll in- 
flation, 0.9 < n(ko) < 1.1. The resulting constraints 
for n(feo) = 1 are displayed in Fig. 2. For n(ko) = 0.9 
and 1.1 the constraints are weakened or strenghtened, 
respectively, at the order of 2 percent. This indicates 
that, for slow-roll inflation models, the constraints are 
not particularly sensitive to the exact shape of the power 
spectrum in the vicinity of the scale of interest. The large 
scale constraints (small k) come from various astrophysi- 
cal sources such as Milky Way disk heating, wide binary 
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FIG. 2: Generalised constraints on the amplitude of the power 
spectrum of primordial curvature perturbations as a function 
of comoving wavenumber (in units of Mpc -1 ). We have as- 
sumed that the power spectrum is scale-invariant over the 
(relatively small) range of scale which contribute to a given 
constraint. Deviations from scale-invariance consistent with 
slow-roll inflation lead to small changes in the constraints (see 
text for further details). 



disruption and a variety of lensing effects. The small 
scale constraints generally arise from the consequences 
of PBH evaporation, in particular on nucleosynthesis and 
the CMB. These evaporation constraints lead to tighter 
constraints on the abundance of PBHs and therefore the 
primordial power spectrum is more tightly constrained 
on these scales. In general the constraints on the ampli- 
tude of the primordial power spectrum span the range 
V-jz < 1CU 2 — 10 _1 with some scale dependence. 



IV. SUMMARY 

We have compiled, and where relevant updated, the 
observational limits on the initial abundance of Pri- 
mordial black holes. We then translated these limits 
into generalised constraints on the power spectrum of 
the primordial curvature perturbation, taking into ac- 
count the full time evolution of the density contrast. 
The constraints on the amplitude of the power spec- 
trum are typically in the range Vk < 1CP 2 — 1CP 1 with 
some scale dependence. This is slightly weaker than the 
Vtz < 1(T 3 - 1CT 2 assumed in Ref. [Hj]. These more 
accurate generalised constraints could be used to more 
accurately constrain the parameter space of slow-roll in- 
flation models (c.f. [H El El EI El H HI ) . 
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APPENDIX A: DENSITY CONTRAST 
CALCULATION 

The gen eral, scalar, perturbed metric can be written 
as (c.f. |8f|): 

ds 2 = a 2 (r){-(l + 2(/))dT 2 + 2B i i s) dTdx l 

+ [(1 - 2^ + -IF.,: dstlx'] , (Al) 

where B^ and E44 can be written as 



B 



E. 



-B. 



k^ kj 



E. 



(A2) 
(A3) 



where (f>,B,ip,E are arbitrary scalar functions which 
describe the perturbations on a Friedmann-Robertson- 
Walker (FRW) background. Performing a first order co- 
ordinate change, f = r + £°, x l = x 1 + £*, the scalar 
metric variables in the new gauge (denoted by a tilde) 
are then given by 



t = <p - e' - ^ q , 
b = B-e+k?, 

E = E-k£, 



(A4) 
(A5) 
(A6) 
(A7) 



where the scalar part of £ l (r, x l ) is defined as = 
— tt-£, primes denote derivatives with respect to confor- 
mal time r, and h = a' /a = aH. The density contrast 
and the velocity perturbation transform as 



5 = 5 + 3/i(l + «;)£ , 
V = V + f, 



(A8) 
(A9) 



where V is related to the scalar part of the 3-velocity 
vector, z/ (s) = ~^V. 

Ref. [88J calculated the evolution of the density and 
velocity perturbations in the conformal Newtonian gauge 
(which has B^ — En = 0). 

They found that during radiation domination (w = 
1/3) for a fluid with vanishing anisotropic stress (and 
hence — ^n), the remaining perturbations evolve 
according to 



4>N 



h («) 



C, 



(A10) 
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s N = (2^ - io(«) - «ji(/c) j r . i \ 1 1 1 



n = (ii(«) - fio(«)) c, 



(A12) 



where k = kj\fiaH and C is a normalisation constant. 

For the PBH abundance calculation we need the den- 
sity perturbation in the total matter gauge (T), which is 
defined by B T + Vr — 0, and E T = 0. Since the comoving 
curvature perturbation, 1Z, is identical to the curvature 
perturbation in the total matter gauge, ipr (see e.g. HH), 
we get using eqs. (|A4HA7p . 



Using the solutions eqs. (|A10|1 and (IA12|) above we get 



K = 



C 



■Jo(k) , 



(A15) 



which reduces in the small k (large scale) limit to C = 
2^/3TZ, allowing us to replace the normalisation constant 
C with the large scale limit of TZ. Equation (|A14[) then 
becomes, using eqs. (|A11|) and (|A12|) . 
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4>N - tVn 



(A13) 



Similarly, the density contrast in the total matter 
gauge (during radiation domination) is 



4v N 

k 



(A14) 



4 / k 



ji(k/VSaH)U. 



(A16) 
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